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I. INTRODUCTION 



Recently, a new type of polynomial systems called exceptional polynomials has attracted 
attention in mathematical physics. Roughly speaking, they are an infinite sequence of 
polynomials which are eigenfunctions of a second-order differential operator but which have 
a non-zero lowest degree m > 0. They are usually generated by an infinite sequence of 
exceptional polynomial subspaces of a fixed codimension m. According to the definition in 
Refs. [1, 2], a /c-dimensional subspace Mk of the jV-dimensional type A monomial subspace 

Ka/- ■ 

M,cVi^\z]^{l,z,...,z^-'), (1.1) 

is called an exceptional polynomial subspace of codimension m — Af — k or an sub- 
space if the linear space T>2{Mk) of second-order linear differential operators which preserve 
Mk is not a subspace of the linear space D2(V^'*) of those which preserve V^'*, namely, 
V2{Mk) ^ r'2(v]/^)- Various mathematical results on the Xi polynomials and Xi subspaces 
are explored in Refs. [1-3]. 

On the other hand, some exactly solvable quantum mechanical potentials whose eigen- 
functions are expressed in terms of X^n polynomials have been reported in the last few 
years [4-6]. All of them were constructed using the technique of supersymmetric (SUSY) 
quantum mechanics [7-10] and have shape invariance [11]. Those findings imply wide appli- 
cability of such mathematical concepts in various physics and mathematical sciences. Thus, 
it would be interesting to see whether there exist other exactly solvable and/or less restrictive 
quasi-exactly solvable models [12, 13] associated with Xm polynomials and X^n subspaces. 

In this respect, the recent research developments have shown that the framework of M- 
fold SUSY [14-17] is remarkably useful for investigating such issues (for a review of A/'-fold 
SUSY see Ref. [18]). This originates from the important fact that jV-fold SUSY is essentially 
equivalent to weak quasi-solvahility (for the precise definitions of the hierarchy of solvability, 
sec, e.g., Refs. [18, 19]) proved in Ref. [15] (for the equivalence in specific models, see also 
references cited therein). Since weak quasi-solvability includes solvability, on the one hand, 
and shape invariance is a sufficient condition for solvability^, on the other hand, shape 
invariance always implies A/^-fold SUSY, 

(A^-fold SUSY) = (Weak quasi-solvability) D (Solvability) D (Shape invariance). 

Hence, all the shape invariant potentials reported in Refs. [4-6] must have M-iold SUSY. 
The ones whose eigenfunctions are expressed in terms of the Xi-Laguerre or Xi-Jacobi 
polynomials are evidently preserve exceptional polynomial subspaces of codimension 1. On 
the other hand, according to Theorem 1.4 in Ref. [2], there is only one exceptional polynomial 
subspace of codimension 1 up to the projective equivalence and its representative can be 
chosen as type B monomial space 

v^^)[.] = <i,. ^-^.^>. (1.2) 



It should be noted that shape invariance is not a sufficient condition for exact solvability: shape invariance, 
as well as solvability, is a local concept while exact solvability is a global one. For the importance of 
recognizing the difference between local and global concepts, see, e.g., Ref. [20]. 
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This monomial space is first considered in Ref. [21] in the context of the classification 
of monomial spaces preserved by second-order linear ordinary differential operators. It 
was shown that any Hamiltonian which preserves the type B monomial space belongs to 
type B A/'-fold SUSY [17, 22]. Thus, we come to the conclusion that all of the models 
whose eigenfunctions are expressed in terms of the Xi-Laguerre or Xi-Jacobi polynomials 
in Refs. [4, 5] belong to type B A/'-fold SUSY up to the projective equivalence. In fact, the 
rational potential V{x) in Ref. [4], Eq. (8), whose eigenfunctions are written in terms of the 
Xi-Laguerre polynomials coincides with one of the type B AA-fold SUSY models Vj^{q) in 
Ref. [22], Eq. (4.2). Explicitly, they are identical to each other up to a multiplicative factor 
V{x) — 2Vj^{x) with the following parameter relations: 

21 + 1 u u{2H-2l-l) 
h^u, ho^——, R^j^ ^ . (1.3) 

Hence, the framework of A/'-fold SUSY would provide a powerful tool for investigating and 
constructing solvable and quasi-solvable quantum systems associated not only with mono- 
mial spaces but also with exceptional polynomial subspaces. On the other hand, it has 
not been reported yet, to the best of our knowledge, any A/'-fold SUSY where a component 
Hamiltonian preserves any subspace of codimension m > 1. 

In this paper, we construct for the first time a family of quasi-solvable and A/'-fold SUSY 
quantum systems where each Hamiltonian preserves an exceptional polynomial subspace of 
codimension 2. We rely on the algorithmic construction developed in Ref. [17]. The resulting 
A/'-fold SUSY systems turn to include as a particular case the rational shape invariant 
potential whose eigenfunctions are expressed in terms of the X2-Laguerre polynomials of 
the second kind in Ref. [5] . Furthermore, we find, in particular, that the two X2 subspaces 
connected by the A/'-fold supercharge would generate the two different kinds of the X2- 
Laguerre polynomials found in the latter reference. 

We organize the paper as follows. In Section II, we introduce a finite-dimensional poly- 
nomial space which turns out to be an exceptional polynomial subspace of codimension 2. 
Then, we present quasi-solvable operators which leave the latter space invariant. With the 
set of the polynomial space and the quasi-solvable operators, we construct in Section HI a 
new type of A/'-fold SUSY systems by applying the algorithm developed in Ref. [17]. We 
present the pair of A/'-fold SUSY Hamiltonians and the A/'-fold supercharges in closed form. 
As a by-product, we obtain another exceptional polynomial subspace of codimension 2 and 
a set of quasi-solvable operators which preserve it. In Section IV, we exhibit a couple of ex- 
amples of quantum mechanical systems which possess the new A/'-fold SUSY. In particular, 
we show that one of them coincides with the rational potential model whose eigenfunctions 
are expressed in terms of the X2-Laguerre polynomials of the second kind recently found in 
Ref. [5]. Finally, we summarize the results and discuss various perspectives of future issues 
in Section V. 



II. AN X2 POLYNOMIAL SUBSPACE AND ITS QUASI-SOLVABLE OPERA- 
TORS 

Our starting point is to consider a vector space of finite dimension A/", 

= «] = <<^i('^; «)>, (2.1) 
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spanned by polynomials (pniz; cc) of degree n + 1 in 2; as follows: 

^n{z; a) = {a + n- 2)^"+^ + 2{a + n-l){a- + (a + n)(a - (2.2) 

where 0, 1) is a parameter. When a = 1, then (pn{z; 1) — {n — l)z"'~^^ by definition. 
Hence, the vector space V)^^"' [z; a] in (2.1) reduces to a space which is equivalent to type 
A monomial space V^2i (cf., Ref. [17]), 

V^^-) 1] = z'{l, z^-') = V]^\ [z]. (2.3) 

Similarly, when a = 0, then (^„(z;0) = (n — 2)z"'^^ — 2{n — l)z" by definition and, in 
particular, 2(fi{z;0) = ip2{z;0) = —2z'^. On the other hand, we can inductively prove the 
following formula: 

r\n—k 

Hence, in the case of a = we have 

V_g'^'')[^; 0] = z^ Vm-1 = z\l, z^-Az,..., z^-^ -{N- 1)2^-2^). (2.5) 

By Proposition 2.5 in Ref. [2] with n = - 1, A = 1, // = 0, and = -j2^-^ the 
{M — l)-dimensional linear space VV-i C V^^f-z] introduced above is a polynomial subspace 
of codimension 1 and its fundamental covariant qv^-ii^) given by 

qvj,_Az) = -W-l){z-2)^-\ (2.6) 

That is, it has a root of multiplicity M — 2 = n — laXz = 2 and thus the space VV-i is an 
exceptional polynomial subspace of codimension 1 which is projectively equivalent to type 
B monomial space [2] 

Vm-. - vi?!iN = (1, . . . , z^-\ z^-'). (2.7) 

Therefore, the constraint a ^ 0,1 prevents the space Vj^j- ^"' [z; a] from reducing to the well- 
studied monomial spaces of types A and B (for the latter issues, see Refs. [1, 2, 21-24]). In 
addition, we can assume Af > 2 without any loss of generality. Indeed, when J\f — 1,2, the 
space V^^^^z; a] is also essentially equivalent to type A monomial space as 

V?'^\z;a] = Mz;a)vi^\z], vt'^\z;a] = Mz;aM^\w], (2.8) 

where w = (p2{z\ a)/(pi{z; a). Thus, we hereafter assume M > 2. 

We shall look for a vector space of linear differential operators of (at most) second order 
which preserve the polynomial space V^. One of the characteristic features of the set of 
polynomials (pn{z;a) defined in (2.2) is the following factorization under the action of the 
operator 9^ — 1: 

(a, - l)(pn{z- a) = -[{a + n-2)z-{n-l){a + n)]z''-^ f {z; a), (2.9) 
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where the common factor f{z; a) is an n-independent polynomial of second degree in z given 

by 

f{z;a) ^ z'^ + 2{a-l)z + {a-l)a. (2.10) 

Taking into account the factorization property (2.9), we easily see that any (at most) second- 
order linear differential operator J of the following form: 

J = P4{z)d^+P3{z)d,+P2{z) + (a, - 1) , (2.11) 

where Pi{z) (i = 1, . . . , 4) are all polynomials in z, maps all the polynomials ipn{z; a) (n = 
1,2,3,...) to other polynomials which are not necessarily elements of V^. The latter fact is 
rather considered as a necessary condition for an operator to preserve the vector space V^. 
Hence, we restrict ourselves to considering linear operators of the form (2.11). As we will 
show shortly, quasi-solvable operators of the form (2.11) which preserve Vj^ = V)^^"' [z; a] C 
"KaA+-2['^] consist of those which have non-trivial pi{z), and thus do not preserve type A 
monomial space V^_|2W' namely, V2{V^^°''^) % 'D2(V^+2)- Hence, the polynomial subspace 
S?'^'^'^\z\ q\ belongs to an exceptional polynomial subspace of codimension 2 which is denoted 
by X2 in Ref. [2]. For > 2, we find that there are four hnearly independent such quasi- 
solvable operators which leave the space Vjj invariant. The first operator denoted by J\ is 
given by 

Ji = zdl - (z - a + 3)a, + + (a, - 1) . (2.12) 

The action of Ji on (pn{z'i a) {n = 1,2,3,.. .) reads as 

Ji<fn{z; a) = -{n + l)(pn{z] a) + (n - l)(a + n)(pn-i{z] a). (2.13) 

Hence, it preserves not only for a specific value of N but also an infinite flag of the 
spaces 

VfcV^-c-.-cV^c---. (2.14) 

Therefore, the operator Ji is not only quasi- solvable but also solvable^. The second operator, 
denoted by J2, which preserves Vjf is given by 

J2 = [z' + (a - 1)(q; + - l)]d^ - \z^ + (A^ + \)z + (a - 1)(3q; + 3A^ - l)\d, 

+ + 1), - 4(„ - 1) + ^ - + <" + ^ - ^' (ft - 1) . (2.15) 

The action of J2 on (pn{z\ a) {n = 1,2,3,.. .) reads as 

(a -|- n — 2)(q; -|- n — 1) J2<^n(^; «) = — (^ — M){a + n — 2)^(^„+i(2;; a) + s~(n, a. A/") 

X (pn{z]a) - (n - l)(a - l){a + M - l)[3a^ + 6(n - l)a + 3n^ - Qn + A\(pri-i{z;a) 
+ {n-l){n-2){a-l){a+M -l){a + n-l){a + n)(pn-2{z-a), (2.16) 



^ This significant characterization of solvability first appeared in Ref. [25] but unfortunately without suffi- 
cient appreciation of the difference between solvability and exact solvability. 
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where s {n, a, Af) is given by 

s- {n, a,X)=2a^ + [n^ - (A/" - 4)n - - 7]a^ 

+ [2n^ - {2U + l)n^ + {U - 12)n + SAT + 9]a 

+ - (TV + 3)n2 + (2A/' - 1)71^ + (A/" + 9)n - 4(A/' + 1). (2.17) 

Hence, it certainly satisfies J2^Jf C Vjj- only for a specific value of A/" G N. The third 
operator, denoted by J3, leaving invariant is given by 

J3 = + 2a + A^ - 1)^2^2 + [{^a-U - 2)z'^ + [3a^ + (X - 2)a - 2{M + l)]z 
+ 4(a - + AT + l)]d, - (AT + - 2)z 

The action of J3 on (^n(-2; ct) {n = 1,2,3,.. .) admits a relatively simple form as 

h(pn{z; a) ={n-Af){a + n- 2)(^„+i(2;; a) + { (3n + l)a^ + \2r? + (AT - 4)n 

+ 3A^ + 2\(x +{Af- l)n{n - 1) - 4(Ar + 1)}^„(2; a), (2.19) 

and thus it also satisfies JsV^ C Vj^ only for a specific value of A/" G N. The fourth 
quasi-solvable operator denoted by J4 has the most complicated form as follows: 

(2q; + A^ - 1) J4 = [{2a + Af-l)z + 3o? + (3A^ - 2)a + A^(A^ - \)\z^dl 

- Bl{z-a,M)d,^M{H ^\)(2a^M -\)z^ + (A/'+ l)(a - l)[3a' 

+ 2(3AA - 8)a + 2(AA - 4)(AA - 1)]^ + 4(a - xf^^^^^^^ {d, - 1) , (2.20) 

where B'^[z]a,N') and D^{z;a,N') are, respectively, given by 

B^(z; a,N) = 2J\f(2a + J\f- l)z^ + [3a^ + (9J\f - 19)a'^ + (5J\f^ - 18A^ + 24)q; 
+ (A^ - 1)(A^2 - 2A^ + 8)1^2 + {a- l)[7a^ + (14A^ - 15)a^ 
+ {3J\f + 1)(3A^ - 8)0; + 2(A^ - 4)(A^2 _ ^^^^^ ^ _ i^^2a^ ^ (5^;^^ ^ ^^^3 

+ {AAf + 3Af+ 13)q;2 + (A^=^ + 2Af'' -7Af- 36)a - 8{Af - l){Af + 2)] (2.21) 

and 

D^{z; a, Af) = [a^ + {2N + 3)a^ + (A^' - - 16)q; - 4(A^ - 1)(A^ + 2)]z 

+ q;[q;=^ + (2A^ + 3)o;' + (A^' - 2A^ - 9)a - 2{M - 1)(A^ + 2)]. (2.22) 

The action of J4 on (pn{z; a) {n = 1,2,3,.. .) reads as 

{a + n — l){a + n)Ji(pn{z'^ a) = {n — M){n — M + l){a + n — 2){a + n — l)(pn+2{,z\ a) 

- {n - M)t];{n,a,J\f)i^n+i{z;a) - (a - l)t^(n, a, A/')(^n(^; a) 

- (n - l)a{a - l){a + - l){a + J^){a + nfifin-i{z; a), (2.23) 
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and thus it surely preserves the space Vj^ for a specific value of jV^ e N. In the above, 
ti{n,a,N') and f^in^a^M) are, respectively, given by 

(2q; + - l)ir (n, a, N) = 2>q^ + (3n + GA/" - 20)a^ - [Sn^ - (9^/" - 26)n 

- 2A^2 ^ 27M - 50]q;3 _ [3^3 + An'^ - {2,U^ - 2,bU + 52)n + TA^^ - AQM + 52]q;2 

- [(STV - 2)n^ + (7A^ - 10)n^ + 2(^2 - 21J\f + 18)n - 2(A^ - l){bU - I2)]a 
-{N-l){n- l)[Un' + 2{U - 2)n - A{Af - 1)] (2.24) 



and 



{2a + Af -l)t^{n,a,Af) = {7n + l)a^ + 2[7n'^ + {7Af - 8)n + 2Af + 6] 



a 



+ [7n^ + 28{Af - l)n^ + {QAf^ - 19Af + 48)n + SAf^ - 19Af - 49] 
+ [(lAAf - ll)n^ + 18(X^ -2U + 2)n^ + {2U^ - 7J\f'^ +X- 97)n 
- IQAf + 7Af + 52]a^ + [{Af - l){9Af - A)n^ + 2{2Af - 7)(A^' + 3)n' 

- {2Af^ + IBAf + ISAf - 70)n - 2{M - - 5Af - 8)]a 

+ 2{Af-l)n{n-l)[Af{Af-l)n-4{Af + 2)]. (2.25) 

Hence, the most general quasi-solvable operator H~ of the form (2.11) which preserves the 
space (jV > 2) is given by 

4 

H~ = - ^ajJj - Co 



i=l 



B{z) + 



4(a - l)Diz) 
f{z;a) 



where (i = 1, . . . , 4) and Cq are constants, while A{z), B{z), C{z), and D{z) are polyno- 
mials in z given by 

"3a;2 + {3Af - 2)a + Af{Af - 1) 



A{z) = a4^z'^ + 



2a + Af-l 
+ aiz + {a — l){a + — l)a2, 

^(^) = - 2«+7^'^j ^^ + [(a - AT - 2)a3 - a^jz 



04 + aa 



z^ + [{2a + Af - l)a3 + a2\z'^ 

(2.27) 



+ {[3a^ + {Af- 2)a - 2(AA+ l)]a3 - {Af + l)a2 - ai}z 

+ 4(q; - 1)(q; + + 1)03 -{a- l){3a + 3J\f- 7)a2 + {a - 3)ai, (2.28) 



C{z) ^A/{A/+l)a4Z^ + {A/+l) 



^,3a^ + 2{3J^-8)a + 2{Af-A){Af-l) 

[a — 1) ; — ^ ; 04 



2q; + A^- 1 



[a - 2)03 + a2 



Z + Co, 



D{z)^ 



{a-l)D^(z;a,J\f) 



tti - [{a^ + Ma -2J\f - 2)a3 + {2a + M - 8)02 - ai]z 



2a + M -I 

- a{a - l){a + N + 1)03 - {a - l){2a + N - l)a2 + aai. 
The operator H becomes solvable only when 

04 = 03 = 02 = 0. 



(2.29) 

(2.30) 
(2.31) 
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III. CONSTRUCTION OF AT-FOLD SUPERSYMMETRIC SYSTEMS 



Now that we have constructed the set of quasi-solvable operators which leave the space 
Vjj- invariant, we are in a position to construct jV-fold supersymmetric systems associated 
with the latter space by applying the systematic algorithm developed in Ref. [17]. The first 
step is to construct an A/'th-order linear differential operator which annihilate the space 
Vjj- and which has the following form (cf., Eq. (2.29) in Ref. [17]): 

PM-^'(ir(^+i:^ri^)^X (3.1) 

\ fe=0 / 

where z{q) denotes the change of variable connecting the variable z with a physical coordinate 
g, and z'{q) is the first derivative of z{q) with respect to q. To construct the operator P^, 
wc first note that the following first-order linear differential operator A" [a) plays a role of 
lowering operator for all the polynomials <^„(2;; a), 

where f{z; a) is the polynomial introduced in (2.10) and the prime denotes derivative with 
respect to z. In fact, the action of A~{a) on ipn{z; a) for an arbitrary n = 1, 2, 3, . . . reads 
as 

i"(a)(^„(2;; a) = {n - l)(^n-i{z; a + 1). (3.3) 
Using the latter formula, we easily prove the following formula by induction: 

f'm-i \ r( ) 

(fn{z;a) — — -(pn-m{z] + m), (3-4) 

J. I 77/ Tfl/ j 



,k=0 



where F denotes the gamma function and the product of operators is defined by 

ki 

YlAk = Ak,Ak,-i...Ak,. (3.5) 

k = kQ 

Then, we immediately know that the properly normalized operator Pj^ of the form (3.1) 
whose kernel is the polynomial space is given by 

'-'(''"/(So 

_ // NAT fiz;a) f{z;a + k + l) f d f'{z;a + k + l) \ 

^"^^ f(z;a+Af)l\ f{z-a + k) \dz f(z;a + k + l)J- 

The next task is to calculate the coefficient w^j!^-^^{z) of the (A/" — l)th-ordcr differential 
operator d^~^ in defined by (3.1). For the latter purpose, we first derive a formula for 
the calculation of w^^i{z) in a general setting, 

n + ^^^^) = n /^(-) di^ + *JJ^i^^)di^ + ■ ■ ■ • (3.7) 

fe=0 ^ ^ \fe=0 / ^ ^ 
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It is an easy task to show inductively with respect to M the following recursion relation: 



Its general solution is given by 



fe=0 



fk{zy 



N-l 



k=0 



gk{z) + {M-l-k) 



f_M 

fk{z) 



(3.8) 



(3.9) 



In our present case, we read from (3.2) and (3.6) that the functions fk{z) and gk{z) are given 

by 



f{z;a + k) 

Then, for the latter fk{z) and gk{z), we have 



f{z;a + k + l) 



(3.10) 



E 

fe=0 



fLiz) _ d 
fk{z) dz 



/M-\ 



/(z;a + fc + l)\ f\z;a^M) f'{z;a) 



,k=0 



f{z;a + k) 



f{z-a+N) f{z;a) 



(3.11) 



and 



fc=0 



fk{z) dz 



^j\f-i 



1- n 



, fc=0 



f{z-a + k + lf 
f{z; a + k)^ 



-(\r . JX^;« + Ar) — /(^;a + A:) 
^ ^/(;.;a+Ar) 2^^f{z;a + k)' 

Substituting (3.10)-(3.12) into (3.9), we obtain 



(3.12) 



f'{z;a) f{z;a+J\/) 
f{z;a) f{z;a + Af)' 



(3.13) 



With the obtained function w^^^^{z), the A/'-fold SUSY pair of gauged Hamiltonians H 
and ^+ are expressed as (cf., Eq. (2.45) in Ref. [17]) 



H^^- A{z 



dz^ 
(1±1) 



\f — 2 

-^A'{z) ± Qiz) 



dz 



-C{z) 



(3.14) 



where Q{z) and C{z) in the present case read as 



C{z) = C{z^ 



2 ■ ■ f{z;a) 

4{a-l)D{z) 



f{z]a) 



(3.15) 
(3.16) 
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It is worth studying another hnear space Vj^ preserved by the partner gauged Hamihonian 
. The latter space is characterized by the kernel of another A/'th-order linear differential 

operator P^, namely, through the relation ker = V^. The operator is obtained from 
in (3.1) by (cf., Eqs. (2.30)-(2.32) in Ref. [17]) 



p^^{-irz'{qr-\Pj^rz'{qr^ 

k=0 



(3.17) 



where the subscript T stands for the transposition in the physical g-space. In the present 
case where Pj^ is given by (3.6), it reads as 



-Af-l 



k=0 



^f{z;a- 1) 



dz f{z;a+Af-k) J f{z;a+Af-k-l) 

-Af-l 



f{z;a) 



}={ ')f{z;a + N-l 



f(z;a) 



y 



(3.18) 



where A^{a) is a first-order linear differential operator defined by 



f{z;a- 1) f d f{z;a) 



+ 



f(z;a-2) \dz f(z;a) 



(3.19) 



The vector space annihilated by P^ is obtained by integrating inductively the differential 
equation P^x(^) = 0- The result is 

f{z; a)f{z; a + N)V^ = vjf ^"^[z; a + TV] = (xi(^; « + AA), . . . , xa/-(^; « + AA)>, (3.20) 
where each Xn{z] a) is a polynomial of degree n + 1 defined by 

Xn{z] a) = (a - n)(a - n + l)z^^^ + 2{a - n - l){a — n -\- l){a - l)z^ 



-\- {a — n — l){a — n){a — l)az 



n-l 



(3.21) 



We find that the obtained linear space in (3.20) is such a space on which the operator 
A"*" (a) introduced in (3.19) acts essentially as a lowering operator. Indeed, we easily derive 
the following formula: 



Xniz; a) 



= (n-l)- 



Xn-i{z;a- 1) 



' f{z;a-l)f{z;a) ^ f{z;a-2)f{z;a-l)' 

With repeated applications of the latter formula, we obtain 



(3.22) 



f m—l 



f}^ ' ) f{z;a + Af -l)f{z;a + Af) 



r{n) 



Xn-m{z; a + J\f - m) 



r{n - m) f{z; a+J\f -m- l)f{z; a + J\f -m) 



. (3.23) 
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Then, combining (3.18), (3.20), and (3.23), we eventually have 



,fe=0 



(Xi(^; g + AA), ■ ■ ■ , xm{z\ a + AA)) 



= 0. 



(3.24) 



As the gauged Hamiltonian if+ preserves the vector space defined by the relation in 
(3.20), namely, C V^, it is evident that the linear operator 



fiz; a)f{z; a + U)H+f{z- a + AT)" V(^; a)-\ 



(3.25) 



preserves the polynomial subspace V^^*"^ [z] a +N] C V'ffl^lz], namely, H+V)^^'\z]a+N] C 

V^y^' [2;; a + M]. On the other hand, it immediately follows from the form of given by 

(3.14)-(3.16) that the operator does not preserve the monomial space V^_^2- -'-^ other 

words, r'2(Vif''^) 2 ^52(Vi^j2)- Hence, the linear space Vj^^'''' spanned by the polynomials 
Xn in (3-21) provides another exceptional polynomial subspace of codimension 2. From 
(3.14) and (3.25), the form of the operator H'^ reads as 



u(A) 



where B'^{z) and C'^{z) are, respectively, given by 

A{z)f{z-a) + 2{a-l)D{z) 



(3.26) 



B+{z) = - (AA - 2)A'{z) - B{z) - 2- 
f\z-a + U) 



f{z;a) 



2A{z) 



f{z;a+J\f) 



(3.27) 



and 



C+(z)=(A^-l) 



'Af-2 



A"{z) + B'{z) 



+ C{z) + 



fiz; a) 



2{2^f-3)A{z) 



+ [{2M - 3)A\z) + B{z)] fiz; a) - 4{a - 1) [D{z) - (A/" - l)D'iz)] } 
2A{z) + [(AA - l)^'(^) + B{z)]f{z; a+Af) 



+ 



f{z;a + U) 
^ A{z)r{z;a)+2{a-l)D{z) 
f{z;a) 



{N-2) 



f{z;a) f'{z;a+Af) ' 
f{z;a) f{z;a + M)^ 



(3.28) 



In the above, the functions A{z), B{z), C{z), and D{z) are given by (2.27)-(2.30). Substi- 
tuting them into the expression (3.27), we obtain 



B+{z)^B+{z;a+NM) + 



f{z;a + J\f) 



(3.29) 
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with 

B+{z; a + M,M) = _ Bt{^:oc+J>fM) ^^ _ ^ ^^^^^ _ ^^^^2 

Za + J\J — 1 

- {[3a^ + {5J\f - 2)a + 2{N - 1)>3 + + 3)a2 - ai}z 

+ 4(q; + A^- 1)(Q; + I)a3 + (Q;+A^- 1)(3Q; + I)a2 -(Q;+A^ + 3)ai, (3.30) 

Dt{z;a+M,X) = , ^ , — -a4- [{a^ + Xa + 2M - 2)a3 

2a + J\J — 1 

+ {2a+Af-3)a2 - ai]z - {a + Af){a + Af - l){a + l)a3 

-{a + J\f-l){2a + Af- l)a2 + {a + J\f)ai, (3.31) 

where 

Bt{z]a + M,M) = 2U{2q^+U ~ l)z^ - [2>a^ - ~ 2{2U'^ + U - 2)q^ 

- 2M^{M - l)]z^ -{a + M-l) [7a3 + 7(AA - l)a' + {2M^ - QAT + 16)a 

- 2{M -1){M- A)\z -{a + M-l) [2a^ + (3jV - 7)q;3 + [M^ - 8J\f - ll)a^ 

- {M^ + 3W - 36)a - 4(AA - - 4)], (3.32) 
£>+ (z; a+HM) = [«' + (A/'+ 3)a;2 + (117V - 16)q; + A{M 2)]z 

+ {a+U)[a^ + (TV + 3)q;2 + (87V - 9)q; + (TV - 1)(37V - 4)]. (3.33) 

Similarly, substituting (2.27)-(2.30) into expression (3.28), we obtain 

C^z)=C-iz;a+M,N) + *-^^P^§^, (3.34) 

}[z]a + M) 

with 

C^{z] a + TV, TV) = M{M + 1)04^' 

'3a3 + 3(Ar - 3)^2 _ (AT^ + 4^/" - 4)a - A/'2(Ar - 1) 



2a+M-l 



-04 



q; 



+ TV)a3 + a2 



^ + 4, (3.35) 



D-(.;a+Ar,AA) = ^ 2. + a|--1 ^^^^^ 

- [{a + N){a + N -l){a + l)az+{a + N -l){2a + N -l)a2 - {a + N)ai]z 

- (a + TV - l){a{a + TV)'a3 + [2q;^ + 3(TV - 1)q; + TV' - TV + 2]a2 
-{a+M)ai], (3.36) 

where Cq is an irrelevant constant which depends on a and TV, and I?^(2;; a+J\f,J\f) is given 
by 

L>+ (z; a + TV,TV) = [0;=^ + (TV + 3)o;' + (8TV - 9)a + (TV - 1)(3TV - 4)]^ 

+ (a + TV) (a + TV - 1) [a' + 3q; + 2(TV - 1)] . (3.37) 
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Prom expressions (3.26)-(3.37), it is evident that the operator consists of four hnearly 

independent operators, each of which preserves the hnear space V^^'^''\z] a + A/"]. Setting aU 
but one of the parameters ai (i = 1, . . . , 4) to be and replacing a + A/" by «• we can extract 
the four second-order linear differential operators which preserve the second exceptional 
polynomial subspace V]^"^^ [z] a\. The first operator Ki associated with oi is given by 

4 

Ki = zd'^ + {z-a- 3)9^ + [{a -l){z + a)d^ + a{z + a - 1)] . (3.38) 

The action of Ki on x„(^; a) {n = 1, 2, 3, . . .) reads as 

KiXn{z] a) = (n + l)xn{z] a) - {n - l){a - n - l)xn-\{z] a). (3.39) 

Hence, it preserves not only V^'^*'\z] a] for a specific value of A/" but also an infinite flag of 
the spaces 

Vf [z; a] C Vf [^; a] C • • • C V^^''^ a] C • • • . (3.40) 

Therefore, the operator Ki is not only quasi-solvable but also solvable. The second operator 
K2 associated with the parameter 02 is given by 

K2 = + (a - 1) (a - A^ - l)]d^ + - (A^ + 2>)z + (a - 1)(3q; - 3A^ + 1)% 
- (AA + l)z - ^Yr--T i -M-2,)z + {a-l){2a-M- l)% 
+ {2a-Af- l)z + 2o? -{Af + 3)a + 2(A^ + 1) }. (3.41) 
The action of K2 on Xniz; a) (n = 1, 2, 3, . . .) reads as 

{a — n + 2){a — n + l){a — n){a — n — l)K2Xn{z; a) = (n — Af){a — n + 2) 

X {a — n + 1)^(0; — n)xn+i{z; a) — {a — n + 2){a — n + l)s~^{n, a,J\f)xn{z; a) 

+ {n - l){a - l){a - Af - l){a - n - l)2[3a^ - 6(n - l)a + 3n^ - 6n + 4]Xn-i{z; a) 

+ {n- l){n - 2){a - l){a - Af - l){a - n)^{a -n- lfxn-2{z] a), (3.42) 

where s'^{n, a,M) is given by 

s+ (n, a, A^) = 2a^ - [n^ - {X - 2)n - X + 3]q;^ 

+ [2n^ - {2J\f + l)n^ - (3A/' + 2)n + A^ + 3]o; 

-n^ + {Af+ l)n^ + {2M + 3)71^ + (AT + l)n - 2{N + 1). (3.43) 

Hence, it certainly satisfies K2V^^'^^\z] a] C V^^'^''\z] a\ only for a specific value of A/" G N. 
The third operator associated with 03 is given by 

^3 = (z + 2q; - A^ - l)z'd^ - {(a + A/')^' + [3q;' - (A^ + 2)a - 2{M - l)]z 

-A{a-l){a-U+ l)]d^ + (A^ + l)az 

- ^1? ~ ^} { [(a' -Na + 2N- 2)z + a{a-l){a-M+ l)]d, 

+ a[{a-Af+l)z + a{a-Af)]}. (3.44) 



13 



The action of K3 on Xn{z; a) admits a relatively simple form as in the case of J3 as follows: 

KsXniz; a) = - {n - J\f){a - n + l)xn+i{z; a) - [(3n + 

- {2r^^+^fn + 3J\f- 2)a + {Af + l)n(n - l)]xniz; a), (3.45) 

and thus it also satisfies K^V^^''\z; a] C V^^''''[z; a] only for a specific value of jV e N. The 
fourth operator K4 associated with the parameter 04 is given by 

(2a-J\f- 1)K4 = [(2a -J\f-l)z + 3a^ - (SX + 2)a + U{U + l)]z^d^ 

- B+{z; a,N)dz + M{N +l){2a-N- l)z^ 

-{Af+ l)[3a^ - 3{2M + 3)a^ + 2{Af^ + 7Af + 2)a - AN{M + 1)]^ 

+ ^7^^ [(« - '^)DUz- a,U)d, + ocDUz- a,U)\ , (3-46) 

J yZ) O!) 

where the functions B^{z; a,J\f) and D^^{z; a,Af) {i = 1, 2) are introduced in (3.32), (3.33), 
and (3.37), respectively. The action of K4 on Xn{z', a) (n = 1, 2, 3, . . .) reads as 

{a - n - 2){a - n - l){a ~ n){a - n + l)K4Xn{z; a) = {n - M){n ~ M + I) 

X (a — 71)^(0; — n + iy'Xn+2{z] a) + (n — N'){a — n + l)^t^(n, Q^,J\f)Xn+i{z\ a) 

-\- {a — l){a — n — 2){q. — n — l)t^(n, «, A/')xn(^; oi) + (n — l)a(a — 1) 

X (a - - 1)(q; - A/')(q; - n - 2)(a - n - 1)^(q; - n)xn-i{z] a), (3.47) 

and thus it surely preserves the space V^^*''* [^; a] for a specific value of jV e N. In the above, 
tf{n,a,J\f) and t^in^a^M) are, respectively, given by 

(2a - AA - l)tt{n, a, AA) = 3q;^ - (3n + GA/" + 20)a^ - [3n^ - {9X + 26)n 

- (2A/'^ + 27A/' + 50)]a^ + [3n^ - ^n^ - {3Af^ + 35Af + 52)n - 7Af^ - 49A/' - 52]a'^ 

- [(3A/' + 2)n^ - {lU + 10)^2 - 2(4Ar2 + 21Ar + 18)n - 2{M + 1)(5A/' + 12)]a 

+ {M+ l)[J\fn^ - (A/" + 4)n^ - 2(3Ar + 4)n - A{M + 1)] (3.48) 

and 

(2q; - A^ - l)4(n, a, A^) = {In + 1)q;^ - 2[7n^ + (7A^ + 3)n + 2X- 11]q;^ 
+ [7n^ + A{lJ\f + 4)71^ + (9A^2 + 7A^ - 42)n + SA^^ - 2lJ\f + l]a^ 

- [(14A^+ll)n=^ + 2(9A^2 + 8A^- 10)n' + (2A^3+A^2 -29A^-7)n 

- (A^ + - 8)]q;' + (A^ + 1) (n - 1) [(9A^ + 4)n' + (4A^' + 7A^ - 10)n 

+ 2A^(A^+l)]Q;-2A^(A^+l)V(n- 1). (3.49) 

It is evident from the resulting actions of i^j {i = 1, . . . , 4) on Xn{z; a) (n = 1,2,3,.. .) 
that the gauged Hamiltonian H^, and thus H'^ as well, are not only quasi-solvable but also 
solvable if and only if the condition (2.31), which is the solvability condition for H~, is 
satisfied. That is, H~ , H^, and get solvable only simultaneously. 

The fact that the four operators Ki {i = 1,...,4) leave the polynomial subspace 
V^y^ invariant in spite of the existence of the fractional coefficients 1/ f{z;a) is par- 

tially explained by factorization properties of the polynomials Xn{z',(y) for all n = 1,2,3,... 
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with the common factor f{z; a) under the actions of two first-order hnear differential oper- 
ators Of and Of, which are analogous to (2.9), 

OtXn{z; a) ^ -{a-n-l){a-n){a-n + l)z''-^f{z; a) (3.50) 

and 

OfXni^', Oi) — [{cy — n){a — n + l)z^ + 2{a — n — !)(« — n + l)(a — l)z 

+ (n - l){a -n-l){a- n){a - l)]z''-^f{z; a), (3.51) 

where Of and Of are given by 

Of = zd, -a, Of = {a- l)d, + z + 2a-2. (3.52) 

It is easy to check that all the fractional parts having the factor 1/ f(z; a) in Ki [i = 1, ... ,4) 
are expressed as linear combinations of Of and Of and thus they map all the polynomials 
Xn{z',o:) to other polynomials in z. The latter fact is inevitable for the operators {i = 
1, . . . , 4) to preserve the polynomial subspace V)y ^ [2;; a]. 

Finally, the jV-fold SUSY pair of Hamiltonians in the physical space is obtained by 
gauge transformations 



z=z{q) 



(3.53) 



where the gauge potentials are given by 



W^{q)^^^\n\2A{z)\± J dz 



(3.54) 



z=z(q) 



2A{z) 

and the change of variable z{q) is determined by 

z'{qr^2A{z{q)). (3.55) 
If we introduce two functions E{q) and W{q) in the physical space as 

following the other types of A^-fold SUSY [17, 22, 24], the A^-fold SUSY Hamiltonians 
are expressed as 

H^^-^-^ + ^W{qr - (^E'{q) - ^^E{qf - 2W'{q) - 2E{q)W{q)^ 

- C{z{q)) ± ^W'{q) + [z"{q)w^^_MQ)) + z' {qfw^^i^:M<D)) , (3-57) 

where C{z) and w^^i{z) are, respectively, given by (3.16) and (3.13). The pair of gauge 
potentials VV^(g) in (3.54) is expressed in terms of E{q) and W{q) as 

>V^(?) = J ^(lE{q) T j dqW{q). (3.58) 
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Similarly, the components of AA-fold supercharges in the g-space are also obtained by 
the same gauge transformations, 



z=z{q) 



Making the repeated use of the identity 



z'{q)-' = z'iq)-''-' — ± F^{q) - kE{q) 



dq 



(3.59) 



(3.60) 



where the function Fa{q) is defined by 



f{z{q)]a) 



(3.61) 



we immediately have 



f{z{q)]a) 



f{z{q)-a+M) 



a+k+l 



(?) + 



Af-1 

n 

fe=0 

M-1 



fiziqy,a + k + l) f d 



f{z{q);a + k) \dq 
^E{q) 



{:r + W(q) 



(3.62) 



and 



-Af-l 



A;=0 



1 - 2k 



E{q) 



X 



f{z{q)-a+M-k) 
f{z{q)-a+N-k-l) 



f{z{q);a) 
f{z{q)-a+Nr 



(3.63) 



It is easy to check that they are connected by the correct relation = ( — l)-^(P^)"'". The 
algorithm automatically ensures [17] that the obtained pair of Hamiltonians in (3.57) 
and the A/'-fold supercharges P^ in (3.62) and (3.63) satisfy the intertwining relations 



PnH- 



H+P, 



H-P, 



(3.64) 



and thus, in particular, the two Hamiltonians H~ and H'^ are almost isospectral. 

By the construction it is evident that the jV-fold SUSY Hamiltonians respectively 
preserve the vector spaces defined by 



-W 



-w 



A/-" 



z=z(q) 



z=z(q) 



{Xi{z{q); a + Af),..., X^f{z{q)■, a + X)) 



f{z{qy,a)f{z{qy,a+M) 



(3.65a) 
(3.65b) 



Hence, if, for instance, both H~ and are Hermitian in a Hilbert space L'^{S) with 5* C M, 
and Vjj- and/or are subspaces of L'^{S), then H~ and/or are not only quasi-solvable 
but also quasi-exactly solvable on S. In the latter cases, the solvable sectors V^^ and/or 
provide parts of the eigenfunctions of H~ and/or if+ defined in L'^{S). 
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A{z) 








Ty|")('s oI' V^(q) 


ai 7^ 0, 02 = 
04 7^ U, oi = 


03 = 
02 = 


03 = 


n 
u 

u 


Rational 


O1O2 7^ 0, 03 


= 04 


= U 




Exponential 


02 7^ 0, oi = 

03 7^ 0, oi = 


as = 
02 = 


O4 = 
O4 = 






Trigonometric or Hyperbolic 


03O4 7^ 0, oi 


= 02 


= 






Other cases 








Elliptic 



TABLE I: The relations between the forms of A[z) and the types of V^{q). 



IV. RESULTING AT-FOLD SUSY PAIRS OF QUASI-SOLVABLE POTENTIALS 

In this section, we shall construct explicitly A/'-fold SUSY pairs of potentials in the 
physical g-space associated with the two exceptional polynomial subspaces V\f^°' and Vv-^*" • 
As in the case of type B A/'-fold SUSY [22] , they have no covariance under the linear fractional 
transformations GL{2, C). Hence, we can at most consider projective equivalence classes as 
was done in Ref. [1] for the Xi subspace, or equivalently, the type B monomial space; the 
forms of potentials are not invariant under all projective transformations in general. In 
this paper, we shall rather content ourselves with exhibiting a couple of particular examples 
since the complete presentation of all the cases would involve many complicated formulas. 
The functional types of potentials such as rational, exponential, hyperbolic, and so on are 
determined by the function A{z) through Eq. (3.55). In Table I, we just show the relations 
between the form of A{z) and the types of the potentials V'^{q). In this paper, we shall 
just consider the two simplest cases where only ai (Example 1) or 02 (Example 2) among 
the four parameters introduced in (2.26) is non-zero. In such restricted analyses, we can 
still reduce the freedom of parameters by considering the scale transformation as was done 
in Ref. [17] for the classification of type C A/'-fold SUSY potentials. That is, by (2.27) and 
(3.55) a rescaling of (i = 1, . . . , 4) and cq by an overall non-zero constant u affects the 
change of variable z{q) as 

z{q; uQi, uco) = z{^/l/q■, a^, Cq). (4.1) 

Prom this and Eqs. (3.54) and (3.56), we easily obtain the scaling relations: 

E{q; utti, z/co) = ^fvE[yfvq\ ai, cq), Fa{q\ uai, ucq) = y/uFaiy/uq; flj, Cq), (4.2a) 

W{q; utti, ucq) = ^/uW{^/uq; ai, Cq), >V^(g; f^Oj, i^Cq) = >y^(0^g; Oi, Co), (4.2b) 

Then, by formula (3.57) the potential terms are scaled as 

V^{q; vai, ucq) = uV^{y/uq; ai, Cq). (4.3) 

Hence, we can fix the value of parameters ai or a2 without any loss of generality. In what 
follows, we shall exhibit for each case the change of variable z = z{q), the functions E{q), 
W{q), and VV^(g) which determine the gauge factors, the pair of potentials V^{q), and the 
pair of solvable sectors V^. They are obtained by the calculations of (3.55), (3.56), (3.54) 
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or (3.58), (3.57), and (3.65), respectively. 



Example 1. A{z) = 2z [oi = 2]. 
Change of variable: z{q) — q^. 
Gauge factors: 

E{q) = -, Wiq) 



2a + AA-8 4(a-l)(g2 + a) 



2q f{q'^]a)q 
>V^(?) = ± ^ In |?| T In |/(g'; a) 



Potentials: 

V-{q) 
V^{q) 

Solvable sectors: 



q^ ^o? - 
~2 ^ 8g2 



+ 4 

2 ■ 8g2 
- a + + 3 - Co. 



1 4(a - l)g2 



+ 



+ 4 



J/ 2 ^2 a + S-Co, 
/(g2;a)2 J 

a-AT+l 4(a+Ar-l)g2 



"Kv = (v5i(5^; «),•••, <^Ar(g^; ot)) 



^a+l/2g-gV2 



= (Xi(g'; « + AA), . . . , XAr(g2; a + AT))- 



^_Q,_A^+l/2gg2/2 

/(g2;a + AA) • 



(4.4) 
(4.5) 

(4.6a) 

(4.6b) 

(4.7a) 
(4.7b) 



Since the discriminant of f{z] a) is 1 — a, both of the potentials V^{q) only have a common 
unique pole at g = for a > 1. Hence, the AA-fold SUSY system is naturally defined 
on the half line q E S = (0,oo). On the latter domain S, it is evident from (4.7) that 
Vj^{S) C L'\S) and Vj^{S) (f. L^S). Therefore, it manifests unbroken M-Md SUSY of the 
system. In addition, the solvability condition (2.31) is satisfied in the system, and thus the 
Hamiltonian H" preserves the infinite flag of the subspaces of L'^{S), 



Vr(5) c V2-(5) c • • • c Vj^iS) c • • • c L'iS). 



(4.8) 



Hence, the Hamiltonian H is not only quasi- solvable but also exactly solvable on S provided 
that the infinite flag constitutes a complete set of the Hilbert space L^{S), namely. 



(4.9) 



The fulfilment of the solvability condition (2.31) also guarantees that the other Hamiltonian 
preserves the infinite flag of the spaces V^^^ {Af — 1,2,3,...) which are not subspaces of 



VnS) C VnS) C • • • C V]^{S) C---(t L'(S). 



(4.10) 



On the other hand, Eqs. (4.6) tell us that the other Hamiltonian has the same form as 
its partner Hamiltonian H". Indeed, if + is identical, up to an additive constant, to H~ with 
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its parameter a replaced by a + M . Hence, the A/'-fold SUSY system has shape invariance 
for all A/" G N. Combining it with the fact that H^V^ C for all U = 1,2,3, .. ., we 
come to the conclusion that both H~ and preserve the two different infinite flags (4.8) 
and (4.10) with suitable parameters. In particular, is also exactly solvable on S if the 
completeness (4.9) is assured, although its jV-fold SUSY sector does not belong to the 
Hilbert space L^(S'). 

Finally, we note that the potential 2V~{q) coincides, up to the scaling factor u = u}/2 
determining the scahng relations (4.1)-(4.3), with case II rational radial oscillator potential 
V^~\x) in Ref. [5], Eq. (2.17), of which the eigenfunctions are expressed in terms of the 
X2-Laguerre polynomials of the second kind L^2n+2{^)- 

Example 2. A{z) = {z^ + C^)/2 [02 = 1/2, = (a - 1) (a + - 1) > 0] . 

Change of variable: z{(i) — Csinhgr. 



Gauge factors: 

£^(g) = tanhq', 

T../ X C ^ 3^ ^ (a-l)(a+A/'-3) 

Wiq) = - cosh q + - tanh q + ^ ^ — ^ 

2 2 ^ cosh 5 

+ [{2a +^f- 3)C sinh q + {a - l){2a + Af - 1)] 



2(a-l] 



Potentials: 



A^- 1±1 



f{( sinh q; a) ( cosh q 



In I coshgl ^ In |/(Csinhg; a)\. 



, e ,2 , AA-l ^ . ^ , 4a^ + 4(A^-4)a+A^^ + 16 

V [q) = — cosh q H — smh q H cq 

8 4 8 



+ 



1 



8 cosh^ q 
2(a-l) 



[4{Af - l)Csinh q + Aa^ + ^{M - 2)a - - 2A^ + 4] 



(sinhg - a - ATH- 3 ^, ^ , 2C sinh g - A/" + 1 

— 2[a — 1) 



/(Csinhg;a) 



/(Csinhg;a)2 



(4.11) 

(4.12) 
(4.13) 



(4.14) 



, , C ^2 3AA-1 . ^ 4a;2 + 4(A^-4)a+A^2^16 

V [q) = — cosh q H Q smh q H cq 

8 4 8 

^ r4(A^ + 1)C sinh q-Aa''- A{Af - 2)a + A^' + 6A^ - 4] 

8 cosh q 



-2{a + Af-l) 
Solvable sectors: 



C sinh q — a + 3 , . . , 2C sinh q + M 
-2ia + N-l) 



/(Csinhg;a+AA) 



f{Csmhq;a+NY 



(4.15) 



(4.16a) 



g-C(smh(/)/2-Cgd(7 

V;^ = (^.(Csinh g; a), . . . , ^^((sinh g; ^)) (,,,hg)A^/^~i/(C sinh g; a) ' 

gC(sinhg)/2+Cgdg 

= (xi(C sinh g; a + AA), . . . , xn{C sinh g; a + -A/")) (,,3hg)A^/V(Csinhg; c + A^) " 

(4.16b) 
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In the above, gdg = arctan(sinh g) is the Gudermann function. The Hamiltonians are 
only quasi-solvable but not solvable since a2 7^ and the solvability condition (2.31) is not 
satisfied. The resulting potentials are of the generalized Poshl-Teller types. 

For a > 1, both of the potentials l^^(q') have no singularities at any finite |g| and thus 
would be naturally defined on the full real axis R. However, neither V^(M) nor V^(M) 
belongs to the Hilbert space L^(]R), which means that both H~ and are only quasi- 
solvable but are not quasi-exactly solvable on M. Hence, j\A-fold SUSY of the system is 
dynamically broken in this example. 

For < 0, that is, for —J\f+1 < a < 1, the change of variable is given by z{q) — \(\ cosh q. 
The form of the potentials in the latter case is similar to the above system for > with 
(, coshg, and sinhg replaced by \(\, sinhg, and coshg, respectively. However, there exists 
a pole at g = in both of the potentials V^{q) irrespective of the value of a. In addition, 
at 2; = 1 — a{> 0) we have /'(I — a;a) — and /(I — a;a) — a — 1(< 0), which means 
that the function f[z; a) has a positive root for all a < 1. Hence, a natural domain of the 
system may be a half line S — (go, 00) with /(|C| cosh go; «) = 0. 

V. DISCUSSION AND SUMMARY 

In this paper, we have constructed a family of quasi-solvable and j'V^-fold SUSY quantum 
systems where each Hamiltonian preserves an exceptional polynomial subspace of codimen- 
sion 2. Wc started with the X2 space Vl;^^° and constructed the four linearly inde- 

pendent second-order differential operators Ji (i = 1,...,4) which preserve it. Wc then 
constructed the A/'-fold SUSY quantum systems by applying the algorithm developed in 
Ref. [17]. As a by-product, we have automatically obtained the other X2 space V^^''''[2;; a] 
and the four linearly independent second-order differential operators Ki {i — 1, . . . ,A) which 
preserve the latter space. This shows one of the advantageous and powerful aspects of the 
framework of A/'-fold SUSY. We presented the two particular examples of the A/'-fold SUSY 
systems. The one is the pair of rational-type potentials which coincide with the rational 
shape invariant potentials in Ref. [5] and thus are not only quasi-solvable but also solvable. 
In addition, it turned out that they admit two linearly independent analytic local solutions. 
The other is the pair of hyperbolic-type potentials both of which are only quasi-solvable. 
Dynamical A/'-fold SUSY breaking would take place in the second example but not in the 
first example. 

The polynomial parts of eigenfunctions of the rational potential V^{q) in (4.6a) are, on the 
one hand, given by the infinite flag of the spaces Vj/^"^ spanned by the polynomials (fn{z; a) 
in (2.2) and are, on the other hand, expressed in terms of the X2-Laguerre polynomials of 
the second kind L^^l^2i^) [^]- Hence, the former polynomials ^n{z', a) and the latter 

^2^u+2i^) iiiust be connected by linear transformations. In fact, the first few X2-Laguerre 
polynomials are expressed by linear combinations of the polynomials <fn{z', a) as follows: 

(«-l)Lg(^)=(^i(^;a), (5.1) 
«-^'2?3 (^) = -'fi2(z; a) + (a + 2)ifii(z; a), (5.2) 
2{a + l)Lg(z) = (fsiz; a) - 2{a + S)(f2{z; a) + {a + 2){a + S)(fi{z; a). (5.3) 

It also indicates that the polynomial system {l'^1+2{^)}i^o would be obtained by the Gram- 
Schmidt orthogonalization of the base system {(^„(^; a)}^^ if it constitutes an orthogonal 
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polynomial system with respect to a certain inner product. 

Similarly, the first few X2-Laguerre polynomials of the first kind L^°j^2(-2) in R-^f. [5] are 
expressed by linear combinations of the polynomials Xni^', (^) in (3-21) as follows: 

a(a + l)LS;^)(z)=Xi(-^;-«), (5.4) 

(a + l){a + 2)Zg(^) =X2(-^; -a) + {a + 3)xi{-z] -a), (5.5) 
2{a + 2){a + 3)L^{;l{z) =X3{-z; -a) + 2{a + 4)x2{-z; -a) 

+ {a + 3){a + 4:)Xii-z;-a). (5.6) 

Hence, we arrive at the following conjecture: 

Conjecture 1 The X2-Laguerre polynomial system of the first kind {-^i"j+2('2)}^o 
would be obtained by the Gram-Schmidt orthogonalization of the base system 
{Xn{—Z', — a)}^i with respect to a certain inner product, while the X2-Laguerre poly- 
nomial system of the second kind {l'^1+2{^)}'u^=o would be obtained by the same 
orthogonalization scheme of the base system {(pn{z; a)}'^^^ with respect to another 
certain inner product. 

Another important remaining issue is to establish a systematic algorithm to calculate the 
characteristic polynomials TV of the superHamiltonian H in the solvable sectors which 
appear in the anti-commutators of the jV-fold supercharges [16], 

{Q^,Q^}=2^^Ar(i?). (5.7) 

In the cases of type A and C A/"- fold SUSYs, it was shown [17, 24] that they are given by the 
generalized Bender-Dunne polynomials of critical degrees [26] and thus are systematically 
calculated via recursion relations. It would be interesting to examine whether a similar 
approach also works for the present AA-fold SUSY systems. 

One of possible continuations of the present work is to construct A^-fold SUSY associated 
with the X2-Jacobi polynomials. In this respect, we note that, in contrast to the X2-Laguerre 
cases where two different kinds were found, the two sets of the Jacobi-type polynomial 
systems associated with the two different extended Scarf I potentials were found to be 
identical with each other [5]. The latter fact together with the present result that the two 
different kinds of the X2-Laguerre polynomial systems would be connected by jV-fold SUSY 
(cf., Eqs. (3.65)) leads to another conjecture: 

Conjecture 2 The pair of solvable sectors of A/'-fold SUSY Hamiltonians associated 
with the A2-Jacobi polynomials -Pit-'+k-^) both spanned by a one common polyno- 
mial system (with possibly different values of parameters) which would generate the 
system {-Pi"+2(-2)}S^o through its Gram-Schmidt orthogonahzation. 

We note that the classification of X2 subspaces has remained unsolved yet. We expect 
that the present work, together with investigation on quasi-solvability and A/'-fold SUSY 
associated with X2-Jacobi polynomials, would provide us crucial clues to it. 

Another possible research direction is to explore quasi-solvability and A/'-fold SUSY as- 
sociated with Xm subspaces for m > 2. It would tell us how the shape invariant potentials 
associated with the Xm {rn > 2) polynomials in Refs. [5, 6] are realized in the more general 
framework of A^-fold SUSY. 
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